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Abstract. We study the resonances of the quantum kicked rotor subjected to an extended initial distribu- 
tion. For the primary resonances we obtain the dispersion relation for the map of this system. We find an 
analytical dependence of the statistical moments on the shape of the initial distribution. For the secondary 
resonances we obtain numerically a similar dependence. This allows us to devise an extended initial con- 
dition which produces an average angular momentum pointing in a preset direction which increases with 
time with a preset ratio. 
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PACS. 37.10.Vz Mechanical effects of light on atoms, 
semiclassical methods. 



molecules and ions - 05.45.Mt Quantum Chaos 



1 Introduction 

The quantum kicked rotor (QKR), can be considered a 
cornerstone in the study of chaos at the quantum level 
PP. Both in theoretical and experimental terms, this topic 
is a matter of permanent attention 2,3,4,5,6,7 . An ex- 
perimental realization of the quantum kicked top (a vari- 
ant of the QKR) has been recently reported [5] and used 
to study the signatures of quantum chaos. The behavior 
of the QKR has two characteristic modalities: dynamical 
localization and ballistic spreading of the variance in reso- 
nance [5] . These behaviors are quite different and have no 
classical analog. They depend on the relationship between 
the characteristic time of the free rotor and the period as- 
sociated to the kick. When the dimensionless period of 
the kick, T, is an irrational multiple of 2ir the average 
energy of the system grows in a diffusive manner for a 
short time and afterwards dynamical localization appears. 
In this case, the angular momentum distribution is char- 
acterized by an exponential decay. When T is a rational 
multiple of 2tt the behavior of the system is resonant with 
a ballistic spreading. In this case the angular momentum 
distribution evolves from the initial one in such away that 
its standard deviation has the time dependence a(t) ~ t. 
This quantum behavior is similar to the optical Talbot 
effect (TUlfTTj in the time dimension. 

The QKR has been experimentally realized through 
a dilute sample of ultra-cold atoms exposed to a one- 
dimensional spatially periodic optical potential that is pulsed 
on periodically in time (to approximate a series of delta 
function kicks) . The quantum resonances and the dynami- 
cal localization have been experimentally observed in refs. 



[T21[T51[T4"lll5] . However with the aim of the present pa- 
per in mind, it is interesting to highlight the work devel- 
oped in refs. [4,5,6,7 , where the influence of an initial 
momentum on the appearance of quantum resonances in 
the QKR is explored experimentally establishing the de- 
pendence of the quantum resonances on the initial velocity 
of the atoms. In ref.[4], the authors showed a sinusoidal 
dependence of the energy on the initial momentum for two 
kicks, and a more complex behavior with the same period 
at four kicks. They also proposed to use the momentum 
dependence of the kicked rotor to select narrow parts of an 
initial momentum distribution. In refs.[5J[7] the dynamics 
of a dilute atomic gas kicked periodically has been stud- 
ied. There, the authors obtain analytical expressions for 
the time evolution of the moments of momentum for ex- 
tended initial conditions, corresponding to the Boltzmann 
distribution of an ideal gas. 

On the other hand, in a recent paper [16], the evo- 
lution of initially extended distributions in the quantum 
walk (QW) on the line was studied. In that work, through 
an analysis of the dispersion relation of the process and 
the introduction of extended initial conditions, continuous 
wave equations are derived. In particular, for a class of ini- 
tial conditions, the evolution is dictated by the Schrodinger 
equation of a free particle. This allows to devise an initially 
extended condition leading to a uniform probability dis- 
tribution whose width increases linearly with time, and 
with increasing homogeneity. Additionally, in some pre- 
vious works [T71H51ll9U20U21[|22[r25] a parallelism between 
the behavior of the QKR and a generalized form of the 
QW was developed, showing that these models have simi- 
lar dynamics. In the same line of thought of reference |16| . 
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in this paper we study the dependence of the QKR in res- 
onance with extended initial conditions. In particular we 
obtain, analytically, a dispersion relation that controls not 
only the velocity of the wavepacket, but also the evolution 
of its shape as a function of time. 

The paper is organized as follows: In the next section 
we present a brief revision of the QKR equations, in the 
third section we present the dispersion relation of the map 
and the group velocity, in section 4 we study the behavior 
of the system under extended initial conditions, in section 
5 numerical results are presented for the secondary reso- 
nances, and in the last section we draw the conclusions. 

2 Quantum kicked rotor 



where the matrix elements of operator U are [Tj 

Ulj =i-U-l) e -ij'eT/h (8) 

Ji being the Ith order cylindrical Bessel function and its 
argument the dimensionless kick strength k = K/%. The 
resonance condition does not depend on k and takes place 
when the frequency of the driving force is commensu- 
rable with the frequencies of the free rotor. Inspection of 
eq. ([SJ shows that the resonant values of the scale param- 
eter r = eT/h are the set of the rational multiples of 2ir, 
t = 2tt p/q. In what follows we assume that the resonance 
condition is satisfied and therefore the evolution operator 
depends on k, p and q. We call a resonance primary when 
p/q is an integer and secondary when it is not. 



The QKR is one of the simplest and best investigated 
quantum models whose classical counterpart displays chaos. 
The Hamiltonian has the following shape 

T 2 00 

U(t) = — +KcoaQ^25{t-nT), (1) 

n=l 

where K is the strength parameter, I is the moment of 
inertia of the rotor, L the angular momentum operator, 
T the kick period and the angular position operator. 
Given the delta function appearing in eq.([T|), the external 
kicks occur at times t — nT with n integer. The index 
n will be equivalent to time in units of T and the two 
operators L and satisfy the canonical commutation rule 

[L,®] = -ih. (2) 

In the angular momentum representation, L|Z) = lh\l), 
the wavevector is 

oo 

\W(n)) = J2 «i(«)IO (3) 

I— — CO 

and the angular momentum and energy are 

oo 

L(n) = (#\L\V)=n E 'l«i(n)| a , ( 4 ) 

/ — — CO 

CO 

E(n) = (V\H\V)=e £ I 2 \ ai (n)\ 2 , (5) 

l— — oo 

where e — h 2 /2I. The evolution operator U is built from 
the Hamiltonian cq.([T]). Then the unitary dynamics of the 
system is described by 

|#(n + l))=U|!P(n)>. (6) 

Substituting eq.([3j in eq.((Sj and projecting over the an- 
gular momentum state |Z), the dynamical equation for the 
amplitudes a; is determined 

oo 

ai(n+l) = E U lj a i( n )' ( 7 ) 

i'=— oo 



3 The dispersion relation and the group 
velocity 

We shall consider the primary resonances p/q — 1 and 
search for the solution of eq.© with the help of discrete 
Fourier analysis. We define the discrete Fourier transform 
of ai(n) as 

S.WE^e'^aifn), (9) 

where uj will be determined as a function of 9. Eq.® can 
be inverted to obtain 

ai{n) = ^ j *Mn)e~ i{n+un) d9. (10) 
Substituting eq. (fTU|) in eq.© and using the identity [2"4"] 

oo 

e - iKCOS e = i-ie-*' Jj(k), (11) 

j=—oo 

the following dispersion relation between to and 9 is ob- 
tained 

U! = K COS 9. (12) 

Therefore, the group velocity of the system is 

Eqs.(|12p and (fTS]) allow us to make simple but relevant 
predictions about the QKR dynamics. In particular, from 
eq. ((T3)) it results that the group velocity will be maximum, 
v g = k, if 9 = £ vt and minimum, v g = 0, if 8 = mir 
with m = 0,±1,±2,.... If the initial state is a wavepacket 
close to some eigenfunction 

Ol (0) = f(l)e ie ° l , (14) 

where f(l) is a smooth function of I and the value of 9q, 
the initial angular position, determines the velocity of the 
wavepacket. If 9q = a sufficiently extended wavepacket 
should stay at rest because v g = 0, while if #o = ±7r/2 it 
should move with maximum velocity v g — ±n. In order 
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to verify this statement, in the next section we choose the 
initial condition as 



oi(0) 



1 P 

ex P (-ttt) 



ctov 2n 



2(7n 



JBqI 



(15) 



where o§ is the standard deviation of the initial distribu- 
tion. We shall then calculate the average momentum and 
energy obtained from this initial condition. 



4 Moments with extended initial conditions 

In the case of the primary resonances, eq.© is solved [T] 
using the recursion relation satisfied by the Bessel func- 
tions. The general solution can be written as 

oo 

ai(n)= J2 (-')'" j ni(^i-iW. (16) 

j=—oo 

where Oj(0) are the initial amplitudes. The probability 
distribution at time n is given by 



Pl(n) = \ai(n)\ 2 , 
that can be expressed as 



(17) 



Pi{n) = J2 {^t 3 a k (Q)a*{Q)J^ k (Kn)J^ 3 {nn). 

j.k— — oo 

(18) 

All the moments of Pi(n) can be calculated analytically 
from eq. ([T8|) using the properties of the Bessel functions. 
We calculate the first and second statistical moments 



Mi(n) = J2 



3 = -oo 



M 2 (n)= fMntf 

j=—oc 



(19) 
(20) 



obtaining: 

oo 

Mi(n) = -ku J2 [o 3 -(0)oj_ 1 (0)] + M x (0), (21) 



M 2 (n) = 



(kti) 2 



l- J2 ^h(QK- 



+ K n (2j + l)3[ aj (0)a; +1 (0)] 

j=—oo 

+M 2 (0), (22) 

where 3? [a;] and 3 [x] are respectively the real and imag- 
inary part of x. We note that the first and the second 



moments are related to the average angular momentum 
and the average energy of the system respectively 



Mr(n) 



M 2 (n) 



L(n) 



E(n) 



(23) 



(24) 



Substituting eq.(JTSJ) into eas. (|21l l2"2"]) and using that 



j = -oo 



J2 ex p(-^2 ) - V2^a, 



(25) 



for a > 1, the following analytical results for the moments 
are obtained 



M\(n) = — nn sm0Q e 8ct o , 



M 2 (n) 



(nn) 2 



cos26> e" 



(26) 



M 2 (0), (27) 



and the variance, a 2 = M 2 — Mf, is then 



(kti) 



2 r 



1 — cos 2#o e~ 



2 snT & 



'o e 



(28) 



Eq. (|27p shows a sinusoidal dependence of the energy with 
0o, the angular initial position for the initial distribution. 
This dependence seems to be confirmed in the experi- 
mental realization of ref.[3]. From the time dependence 
of eq. (|28p , it is clear that the ballistic behavior of the res- 
onance is independent of the initial condition. However, 
from this equation we also conclude that, for n fixed, 
a (n) i— !• (To 2 when <jq 3> 1. This means that when uq 
is very large the system can take a large time to exhibit 
its ballistic behavior. Therefore, the experimental observa- 
tion of this behavior will be subordinated to the extension 
of the initial distribution. This result could explain the 
absence of ballistic growth in some experimental realiza- 
tions [13] . Additionally, our previous reasoning can also be 
related to refs.[51[7] where the authors find analytical ex- 
pressions for the moments with extended initial conditions 
that correspond to the Maxwell-Boltzmann ideal gas dis- 
tribution at finite temperatures. The parameter oo m our 
formulation can be qualitatively related to the tempera- 
ture, a high temperature corresponding to a high value of 
<To and in both formulations the behavior of the moments 
is consistent. 

From ea. (|2"6"]) and ea.([2"5]) the QKR group velocity is 
calculated 



dMx{n) _ 1 dL{n) 



dn 



h dn 



— — /tsmf/o e 



8"0 



(29) 



and we confirm that if 9q = then v g = 0, while if 
6q = ±7r/2 then v g ~ ±k for do > 1. This result is very in- 
teresting because it shows that if one chooses the value of 
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Fig. 1. Using eq.flU) the function 2- 



is presented, in a 
log-scale, as a function of the initial angular position, 9o/ir, for 
three values of ao'- 1, 10 and 100. The calculation corresponds 
to the primary resonance p/q = 1. 



#o appropriately in the initial condition of the QKR, one 
can get an average angular momentum pointing in a pre- 
set direction and increasing in time with a preset ratio v g . 
In order to visualize the temporal evolution of the prob- 
ability distribution in momentum space, we present the 
figH] obtained numerically using the original map eq.©. 
There, the initial conditions 6q — tt/2 in ea.([T5|) have been 
chosen so that the group velocity is negative, hence, the 
first moment moves to the left. The spreading of the wave 
function in the angular momentum space as time elapses 
can also be appreciated in the figure. 

Strictly, from eg. ([25]) and fig. [1] the values of 0q and 
(To cannot be chosen in order to produce a{n) = ao with 
v g 7^ 0. However, for 9q = ±ir/2 and ao > 1 it is possi- 
ble to produce a wavepacket with a constant shape that 
travels with constant velocity for a very long time. Addi- 
tionally, in the asymptotic limit ao — > oo, from eci. ([2l))) we 
recover the expected equation eq. (fK3"| that corresponds to 
the group velocity for a train of plane waves. 

Therefore, we have demonstrated that the dispersion 
relation eq. (fT3"]) controls the velocity of the QKR wavepacket 
and the temporal evolution of its shape, see ecis. (j26l l27|) . 



5 Secondary resonances 

We now study the incidence of extended initial distribu- 
tions on the secondary resonances. In this case p/q is not 
an integer and eq.© becomes too complicated to han- 
dle analytically, an exception being the antiresonance case 
p/q = 1/2 [S]. Eq.© is no longer invariant under transla- 
tions in the angular momentum. This invariance was fun- 
damental for the previous analytical treatment. We there- 
fore study numerically the influence of the extended initial 
conditions in the dynamics of the secondary resonances, 




Fig. 2. The probability distribution in momentum space as a 
function of the dimensionless angular momentum is presented 
for three times t = 0, t = 250 and t = 500. The figure was 
obtained numerically from the map eq.© with k = 0.25 in the 
primary resonance p/q = 1. The initial conditions have been 
given by eq. (|15|) with ao = 1 and 9o = n/2. 




Fig. 3. The first statistical moment as a function of the di- 
mensionless time n for two values of n = 0.25 and 0.4. The 
dashed lines are lineal adjustment, the angular initial position 
is 8o = 7r/2 and ao = 1. The calculation corresponds to the 
secondary resonance p/q = 1/3. 



using eq.© with p/q = 1/3, and initial conditions given 
by ea. lfTS")) . 

In fig. [3] we present the numerical calculation of the 
first statistical moment Mi(n) as a function of time n, 
for two different values of n with 8q = tt/2 and ao = 1. 
This figure shows a complex oscillatory behavior where it 
is possible to define an average lineal behavior, shown by 
the blue dashed lines. The slopes of these straight lines de- 
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fine the corresponding group velocities for the secondary 
resonances. In fig. [4] we present the numerical calculation 
of the group velocity of the system as a function of the an- 
gular initial position. This figure shows a quasi sinusoidal 
dependence of v g with 0q, that reminds us of the behavior 
that was found for the primary resonances, eq. (p?5|) . How- 
ever, if this behavior is compared to that of the previous 
case, now the oscillation frequency duplicates its value and 
its amplitude decreases by some orders of magnitude (note 
that this curve in fig. [4] is multiplied by 60). Addition- 
ally, fig. [3] also shows that v g has a dependence with the 
strength parameter k. We have obtained numerically this 
dependence and we synthesize the result in fig. [5] From 
this figure, it is clear that the simple proportionality ob- 
tained for the primary resonance, see ea.(|26p. has been 
lost. 




-1 -0.5 0.5 1 

e,A 

Fig. 4. The dimensionless v g as a function of the initial angular 
position 9q is presented. A thin red line is used for the primary 
resonance and a thick black line for the secondary resonance 
p/q = 1/3. Both calculations have <tq = 1 and k — 0.25 but v g 
for the secondary resonance is multiplied by a factor of 60. 




0.5 1 



K 

Fig. 5. The dimensionless v g as a function of the strength 
parameter k for the secondary resonance p/q — 1/3 with cto = 1 
and 8 = tt/4. 

Additionally, from the analytical treatment we show the 
sinusoidal dependence between the energy and the initial 
average position of the system. 

We also study, numerically, the influence of an ex- 
tended initial distribution in the behavior of the secondary 
resonances of the QKR. We obtain a relation between the 
parameters of the extended initial condition and the dy- 
namical average of the system. We show, as in the case of 
the primary resonances, that if one chooses the value of 
9q appropriately in the initial condition, one can get an 
average angular momentum pointing in a preset direction 
and increasing in time with a preset ratio v g . 

We acknowledge the support from PEDECIBA, CSIC, 
ANII and thank V. Micenmacher, E. Roldan and G. J. 
Valcarcel for their comments and stimulating discussions. 
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6 Conclusion 

In summary, through the use of the discrete Fourier trans- 
form in the study of the primary resonances of the QKR, 
it is straightforward to get an explicit dispersion relation. 
This dispersion relation is a powerful tool for predicting 
the QKR dynamics because it controls not only the veloc- 
ity of the wavepacket, but also the evolution of its shape 
as time elapses. We show that the extension of the ini- 
tial distribution can delay the characteristic time that the 
resonant evolution needs in order to exhibit its ballistic be- 
havior. This result could explain the absence of the ballis- 
tic growth in some experimental realizations of the QKR. 
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